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Abstract. We consider the fractional Laplacian on a domain and investigate the asymp- 
totic behavior of its eigenvalues. Extending methods from semi-classical analysis we are 
able to prove a two-term formula for the sum of eigenvalues with the leading (Weyl) term 
given by the volume and the subleading term by the surface area. Our result is valid under 
very weak assumptions on the regularity of the boundary. 



1. Introduction and main result 

1.1. Introduction. In this paper we study the asymptotic behavior of eigenvalues for frac- 
tional powers of the Laplacian. The operator (— A)^ with < s < 1 appears in numerous 
fields of mathematical physics, mathematical biology and mathematical finance and has at- 
tracted a lot of attention recently. The key difference between this operator and the usual 
Laplacian is the non- locality of (—A)'', which allows one to model long-range interactions in 
applications and leads to challenging mathematical problems. 

From a probabilistic point of view, the fractional Laplacian of order s on a domain C M'^ 
can be defined as the generator of the 2s-stable process killed upon exiting fi. A more 
operator theoretic definition, which we employ here, is in terms of the quadratic form 

Cs,J' [ ^-4^^-^^dxdy= [ \pnn{p)\'dp, (1.1) 
\x — y\ ^ jRd 

restricted to functions u £ H^(R'^) which satisfy n = in IR'^\0. Here H^(R'^) is the Sobolev 
space of order s, u{p) = {2tt)~^^'^ J e~^^'^u{x) dx is the Fourier transform of u and Cs^ is an 
explicit constant given in (jl.Sp . The identity in (jl.ip is an easy consequence of Plancherel's 
theorem. 

For bounded domains the spectrum of the fractional Laplacian is discrete and we denote 
its eigenvalues (in increasing order, repeated according to multiplicities) by A^'^^ Our main 
result in this paper is a two-term asymptotic expansion of the sum of these eigenvalues, 

j- A(f) = C^^IM~^''^ N^'/'' + cg|ai7pr('^-i+2^)/'^ iV(2«-i)/'i (1 + o(l)) as iV ^ oo . 

n=l 

(1.2) 

Here and \dVl\ denote the d-dimensional measure of 17 and the (d— l)-dimensional surface 
measure of dO., respectively, and C^j and C^^j are positive, universal constants, depending 
only on d and s, for which we shall obtain explicit expressions. Our result is valid for 
non-smooth domains, requiring only that dO. G C^'" for some (arbitrarily small) a > 0. 
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It is remarkable that, despite the fact that we are deahng with a non-local operator, both 
coefficients in (jl.2p have a local form, depending only on 0, and dO,, just like in the case of 
the Laplacian. This will become clearer from the reformulation given in Theorem [T] below. 

In order to avoid confusion, we emphasize that the fractional Laplacian of order s on 
a domain 17 is different from the Dirichlet Laplacian on Q raised to the s-th power. For 
the Dirichlet Laplacian, and hence for its fractional powers, asymptotics analogous to (II. 2p 
are well-known. One of our results is that, while the first terms in (II. 2p coincide for both 
operators, the second terms do not. This means, in particular, that our result cannot be 
obtained from the study of the (local) Dirichlet Laplacian, and that our analysis needs to 
take into account the non-locality inherent in (jl.2p . For further results about the relation 
between the fractional Laplacain on a domain and the fractional power of the Dirichlet 
Laplacian we refer to |CS05j : see also Section [6] below. 

The one-term asymptotics Aj^^ = ^ ci]l\^\~^'^'^ Ar2s/rf(i + ^(i))^ which is a frac- 
tional version of Weyl's law, is a classical result of Blumenthal and Getoor |BG59| . More 
recently, Bahuelos and Kulczycki |BK08j and Bahuelos, Kulczycki and Siudeja |BKS09j 
have shown a two-term asymptotic formula for XlJ^Li 6^P(~*^" ^) as t — )• 0. Note that 
exp(— iAn ■*) and ;^ X]^=i ^ correspond to the Abel and Cesaro summation of the 

(s) 

sequence X-h , respectively. As is well-known, asymptotics of Cesaro means imply asymp- 
totics of Abel means, but not vice versa. Hence for C^'" domains we recover and improve 
upon the result of |BK08llBKS09] . 

This is, actually, a significant improvement since our asymptotics are no longer derived for 
the infinitely smooth function e~''^ of the fractional Laplacian, but, as we shall see shortly, 
for the Lipschitz function (A — £')+. Moreover, since we are no longer able to apply the 
probabilistic machinery available for the partition function, we have to find new and more 
robust tools. Our methods also work for the ordinary Dirichlet Laplacian on a bounded 
domain, and in |FG11| we use the techniques developed here to give an elementary and 
short proof of two-term asymptotics in that case. 

Another point in which we go beyond [BK08llBKS09] is that we give an expression for the 
constant j in (jl.2p in terms of a model operator on a half-line instead of a model operator 
on a half-space. In this way our expression is similar to familiar two-term formulas in semi- 
classical analysis; see, for instance, |SV96) . This is possible due to some recent beautiful 
results of Kwasnicki [KwalOa] about a general class of half-line operators. 

We find it convenient to prove ()1.2|) in an equivalent form, namely 

oo 

^(a-AW) =4;]|17|A1+'^/2--l2|9J7|A1+('^-i)/2^(1 + o(1)) as A ^ oo . (L3) 

n=l ^ 

Here x+ := max{x,0} denotes the positive part of a number x. (The fact that (jl.2p and 
()1.3p are equivalent is well-known to experts in the field, but we include a short proof in the 
appendix for the sake of completeness, see Lemma [20l) Note also that (jl.3p can be rewritten 
as 

oo 

^ (l - /i2^A(f)) =L^^^l\n\h''^-Lf^^\dn\h-'^+\l + o{l)) as/i^O+, (L4) 

n=l 
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and this is the form in which we shall state and prove our main theorem. The small parameter 
h has the interpretation of Planck's constant and (jl.4p emphasizes the semi-classical nature 
of the problem. 

Our approach extends the multiscale analysis to the fractional setting. By this we mean 
that we localize simultaneously on different length scales according to the distance from the 
boundary. Of course, a main difficulty when dealing with our non-local operator comes from 
the treatment of the localization error. At this point we have to improve upon previous 
results from [LY881[SSS10| . Another major impass, as compared to the local case, is the 
analysis of a one-dimensional model operator for which an (almost) explicit diagonalization 
is far from trivial. This is where Kwasnicki's work |KwalOa] enters. It requires, however, 
still substantial work to bring these results into a form which is useful for us. We will explain 
the strategy of our proof in more detail in Subsection 11.31 after a precise statement of our 
main result. 

Throughout this paper we assume that the dimension d> 2. In the one-dimensional case 
(the fractional Laplacian on an interval) considerably stronger results are known [KKMSlOl 
IKwalOb] . The powerful methods developed there are, however, intrinsically one-dimensional 
and seem of little help in the multi-dimensional case. The question raised in |BKS09j of 
whether an analogue of Ivrii's two-term asymptotics |Ivr80| holds for An ^ in d > 2 without 
Abel or Cesaro averaging remains a challenging open problem. 

1.2. Main Result. Let Q, C M'^, d > 2, be a bounded open set. For h > and < s < 1 
let 

Hn = i-h^AY - 1 
be the self-adjoint operator in L'^{Q) generated by the quadratic form 



{u,Hnu) = / {\hp\^' - I) \u{p)\^ dp 

with form domain 

n'in) = |n G H'(M.'^) : u = on M*^ \ . 

For < s < 1 we have the representation 

f f \u{x) — u(y)p f 

{u,Hqu) = Csdh^" / / dxdy- / \u{x)\'^ dx 

jRdJ^d \x-y\'^+'''' Jq 

with constant 

Cs,d = 2 TT ' \T{-s)\ ^ 

Our main results hold without any global geometric conditions on 0. We only require 
weak smoothness conditions on the boundary - namely that the boundary belongs to the class 
C^'" for some a > 0. That is, the local charts of dCl are differentiable and the derivatives 
are Holder continuous with exponent a. 

Theorem 1. Let < s < 1 and assume that the boundary of satisfies d^l £ C^'" with 
some < a < 1. Then 

Tr{Hn)- = h-'' - if} \dn\ h-''+^ + (1.6) 
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with Rh = o{h^'^+^) as h^0+. Here 



(2^) 

(2) 

and the positive constant L^^ is given in h3. 



4;] = 7i/ {\p?'-i).dp (1.7) 



More precisely, we have the lower hound Rh > —Ch '^+1+'=- for any 

< e_ < 



7^2 ifl/2<5<l, 
if < s < 1/2 , 



a+l+2s 

and the upper hound Rh < Ch^'^^^^''+ for any 

a 

< e+ < if 1 - d/4 < s < 1 , 

a + 2 

a(2s-l + d/2) 

< e+ < ^ ^ '.J if < s < 1 - d/4 . 

a + 2s + (i/2 ' 

We do not claim that our remainder estimates are sharp. They show, however, that our 
methods are rather explicit and they correctly reflect the intuitive fact that the estimate 
worsens as the boundary gets rougher. We also mention that for not too small s we (almost) 
get the same remainder estimate /;,-'^+i+°/{"+2) our method yields in the local case 
s = 1 [FGTT] . 

In Section [U] we will derive several representations of the constant ^ in (|1.6|) . One of 
these, which emphasizes the semi-classical nature of the problem, leads to a rewriting of 
([LH]) as 

. /X.^ (N- 1^ - /X.^„ c(i,'r^-,|^ . , ,i.s, 

where T*VL = VLy^W^ and T*dVL = x R'^-i are the cotangent bundles over J7 and dO,, 
respectively, and where da is the surface element of dQ. Here C is a universal (i.e., depending 
on s, but independent of O or d) function, which has the interpretation of an energy shift 
(the integral of a spectral shift). It is given in terms of a one-dimensional model operator 
A'^ on the half-line and its analogue A on the whole line (see Section [3]) by 

/■oo 

C(/i) = /U^^ / [a{t,t, n) — a^{t,t, fj,)) dt , /x > , 

where a{t,u,fi) and a'^{t,u,fi) denote the integral kernels of {A — /i)_ and (A"*" — fj,)^, 
respectively. Another representation, derived in Remark[2l shows that our result is consistent 
with the result of |BK08[lBKS09] . 
In Section [6] we also prove that 

Moreover, we compare this constant with the one obtained from the corresponding fractional 
power of the Dirichlet Laplacian. 

Proposition 2. Let < s < 1 and assume that the houndary of i7 satisfies dil £ C^'°' with 
some < a < 1. Let —Aq be the Dirichlet Laplacian on 17. Then 

Tt{{-h^Any = Lf}\n\h~'' - Lf^Jdnih-''-^' + Rh (1.9) 
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with Rfi = o{h^'^^^) as h ^ 0+. Here L^^^ is the same as in ()1.7p and L^g \ satisfies 

4']<^S- (1-10) 

In other words, the operators Hq^ and (— /i^Afj) — 1 differ semi-classically to first sub- 
leading order. 

1.3. Strategy of the proof. The proof of Theorem [1] is divided into three main steps: 
First, we localize the operator Hq^ into balls, whose size varies depending on the distance 
to the complement of VL. Then we can analyze separately the semiclassical limit in the bulk 
and at the boundary. 

The key idea is to choose the localization depending on the distance to the complement 
of rj, see |Hor85t Theorem 17.1.3] and |SS03j . Let d{u) = inf{|x — u\ : x ^ Q} denote the 
distance of n G to the complement of Q. We set 

l{u) = l(l + [diuf + ll)-'^'Y\ (1.11) 



2 

where 0<Zo^l/2isa small parameter depending only on h. Indeed, we will finally choose 
^0 proportional to h^ with suitable < /? < 1. 

In Section [5] we construct real- valued functions (pu £ C^(M'^) with support in the ball 
Bu = {x : \x — u\ < l{u)}. For all u G M*^ these functions satisfy 

||0«IL<C, ||V<A„|L<C/(n)-i (1.12) 

and for all x G M"* 

(pl{x)l{u)-'^du = 1. (1.13) 



Here and in the following the letter C denotes various positive constants that are independent 
of u, Iq and h. 

Proposition 3. There is a constant C > depending only on s and d such that for all 
< Iq < 1/2 and all < h < C~^lo the estimates 

< TT{Hn)-- [ TricPuHnM^liuyUu < C h-^^"" l^^ R.UKh) 
hold with a remainder 

^1 if 1 - (i/4 < s < 1 



Rloc{h, Iq 



|ln(/o//i)|^/^ if0<s = l-d/4 . 
{lo/hf-^'-'^/^ if < s < 1 - (i/4 



In view of this result, one can analyze the local asymptotics, i.e., the asymptotic behavior 
of TT{(j)uH^(j)u)-, separately on different parts of 0. First, we consider the bulk, where the 
influence of the boundary is not felt. 

Proposition 4. Assume that (p G Cq{Q) is supported in a ball of radius I > and that 

||V(/>||oo < cr^ (1.14) 

Then for all h > the estimates 

-C/°'-2/i-^+2 < Tr{cl)Hn(l))_-L';l [ <i?{x)dxh-'^ < 

' Jn 
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hold with a constant depending only on the constant in ^1.14\ )- 

Close to the boundary of more precisely, if the support of cj) intersects the boundary, a 
boundary term of the order h~'^~^^ appears. 

Proposition 5. Assume that (p G Co(I^'^) ^-^ supported in a ball of radius < I < 1 in- 
tersecting the boundary of $7 and assume that is satisfied. Then for all h > the 
estimates 

-Rhd{l,h) < Tr{(t>Hn4>)--L^^l I ^\x)dxh-'' + Lfl [ <t)\x)da{x)h-''+^ < i?bd(/,/i) 

' JQ ' Jdn 

hold. Here da denotes the {d — 1)- dimensional volume element of and the remainder 
terms satisfy for any < 5i < 1 and < 62 < min{l, 2s} 

/ j^d-l-Si j2a+d-l 
Rhdil, h) < Cs,,52 ^ f^d-l-5r + + + -f^J ' 

with constants depending on 61, 62, ^, \\<P\\oo <ind the constant in |j. 

Based on these propositions we can complete the proof of Theorem [TJ 

Proof of Theorem\^ In order to apply Proposition [5] to the operators cpuHficpu, we need to 
estimate l{u) uniformly. Let 

U{n) = {ueR'^ : BuHdn^^} 

be a small neighborhood of the boundary. For u G [7(17) we have d{u) < l{u), which by the 
definition of l{u) implies 

l{u) < lo/Vs. (1.15) 

In view of (jl.l2p and ()1.15p we can apply Proposition |4] and Proposition [5] to all functions 
^u, -u G M'^, if /q is sufficiently small. Combining these results with Proposition [3] we get 

-C [ l{uy^duh-'^+^ - [ R^A{l{u),h)l{uY'^du 
Jn\u{n) Ju{n) 

du , _^ ^c?:) f f ,0 du ^-d+i 



< TV(i/oL-L« / f <fiSx)dx^^h-'^ + Lf\f [ <ti{x)da{x 

< [ RbA{l{u),h)l{u)-^du + Ch-'^+'^l^^Rxoc{lo,h) . 
Ju{n) 

Now we change the order of integration and in view of ()1.13p we obtain 

-C [ l{u)-^duh-'^-^^ - [ RMiKu),h)l{u)-'^du 
< Tr {Hn)_ - L^] \n\ /i"^ + if] \dn\ h'^^^ 



l{u 



< / RMiliu),h)l{u)-'^du + Ch-'^+^lQ^Rioc{lo,h). (1.16) 
Ju{n) 



It remains to estimate the error terms. 
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By definition of l{u) we have 

l{u)> ^inm{d{u),l) and /(u) > ^ (1.17) 
for all u G M'^. For u £ Q \ U{Q,), we find d{u) > l{u) > /o/4. Hence, we can estimate 

[ l{u)-^du < C ll+ [ d(u)-^du] <C [1+ [ t-^\dnt\dt] , 

Jn\u{n) \ J{d{u)>io/'i} J \ Jio/i ) 

where \drtt\ denotes the surface area of the boundary of fit = {x G $7 : d{x) > t}. Using 
the fact that \dQt\ is uniformly bounded and that \dQt\ = for large t, we get 

[ l{u)-^du < CIq^ . (1.18) 
Jn\u{n) 

For u E U{Q) the inequalities ()1.15p and (|1.17p show that l{u) is proportional to Iq. Since 
BuDdn^^ we find d{u) < l{u) < CIq and 

/ l{u)''du < Cl^ [ du < Cl^+^ , (1.19) 

Ju{n) J{d{u)<io} 

for any a € M. 

We insert (fTTHD and (frT9]l into (fTTBD and get (using the fact that h < C^^/q) 
-C (l-^^h^^ + + q+^h-^) < h''-^ (Tr {H^)_ - \n\ h-'' + L^} \dn\ /i-'^+i) 

< C (V^i/i^i + /o"i/ii?ioc(/o,/i)) • (1.20) 



In order to choose /q we need to distinguish several cases. For the lower bound we recall 
that < 62 < min{l, 2s}. The stated lower bound on follows with Iq proportional to , 
where /3 = (1 + 62)/{l + a + 62). 

For the upper bound we have < 5i < 1. If 1 — d/4: < s < 1, we pick Iq proportional 
to /i^, where /3 = (1 + 61) /{I + a + 61). If < s < 1 - d/4, we pick /i'^, where (3 = 
{2s + d/2)/{a + 2s + d/2). This completes the proof of Theorem [TJ □ 

The remainder of the text is structured as follows. First we analyze the local asymptotics 
in the bulk and prove Proposition [H This is done in Section [2j In Section [3] we consider 
the local asymptotics in the case where is replaced by a half-space. We reduce the prob- 
lem close to the boundary to the analysis of a one-dimensional model operator given on a 
half-line and give an analogue of Proposition [5] for a half-space. In Section H] we show how 
Proposition [5] follows from the previous considerations by local straightening of the bound- 
ary. In Section O we perform the localization and, in particular, prove Proposition [3J In 
the appendix we provide some technical results about the one-dimensional model operator 
introduced in Section [3l 

Notation. We define the positive and negative parts of a real number x by x± =max{0, ±x}. 
We use a similar notation for the heavy side function, namely, = 1 if ±x > and = 
if zizx < 0. For a self-adjoint operator X, the operators X± and are defined similarly via 
the spectral theorem. 
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2. Local asymptotics in the bulk 

This section is a warm-up dealing with the spectral asymptotics in the boundaryless case. 
Although the estimates in this case are essentially known, we include a proof for the sake 
of completeness and in order to introduce the methods that will be important later on. We 
divide the proof of Proposition H] into two subsections containing the lower and the upper 
bound, respectively. The operator 

Ho = {-h^AY - 1 in l2(M'^) , 

defined with form domain H^{M.'^), will appear frequently. 

2.1. Lower bound. The lower bound is given by a variant of the Berezin-Lieb-Li-Yau 
inequality, see |Ber72[lLle73llLY83j . For later purposes we record this as 

Lemma 6. For any cj) G L^(M'^) and h > 

Proof. We apply the variational principle for the sum of the eigenvalues 

-TT{ct)Hn(t>)_ = inf Tr(7</./?n0) , 

0<7<1 

where the infimum is taken over all trial density matrices, i.e., over all trace-class operators 
< 7 < 1 with range belongig to the form domain of Hq^. We apply this twice and find 

Tr(0Ff,,^)_ < T^{(t>Ho<t))_ < Ty:{<P{Ho)_4>) . 

Applying the Fourier transform to diagonalize the operator (Hq)- yields the bound 



T¥(0(Fo)_0) = jj </.(x)^(H^^-l)_ dpdx = L\\>^ J m'dxh 

as claimed. □ 

2.2. Upper bound. We now assume that (j) satisfies the conditions of Proposition HI In 
particular, we assume that (j) has support in VL. To derive the upper bound we put 7 = (ffo)-j 
i.e., 

7(x,y) = {2T^h)-'^ [ e'^^^'-y'^'^dp, 
J\rA<i 

and obtain that 

-Tv{<t>Hn<fy_ < TT{-f<pHn(l>) = Tr(7</./7o0) 

_^(||M^Ar/Ve^W/.||i- 11011^)^. (2.1) 

Lemma 7. For cj) € Co°(M'^) and h > we have 

||(-/j2A)^/20e-W/^||2 = 11,^11^ + I (^^{\p + hv\^' + \p-hr^\^')-\p\' 
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Proof. By Plancherel's theorem we get 



s/2,Jp-/hu2 



dxdyd^ 
{2'KhY 



dxdyd^ 



{27rhr ■ ^^-^^ 

In the first two terms we perform the ^ integration and either the x or the y integration to 
arrive at 



dxdyd(, 



\P\ 



1 2s 



x) dx . 



(2.3) 



We are left with calculating the third term in (|2.2p . Again, by Plancherel's theorem we see 
that it equals 



2s 



We can write 



l-e 



-iz-ri/h 



1 



'■r)/h 



J{p-i)-z/h 



drjdzdS, 



2-e 



iz-ri/h 



and perform the integration in z and ^ to obtain 



2s 



(x) -(/)(y)|2e^(P-«)-(^-S')/^ 



-iz-rj/h 



dxdydS^ 
{2-KhY 



Hence, combining (j2.2p . (|2.3p and (j2.4p yields the claim. 
In view of identity (j2.ip and Lemma [7] we conclude 



di] . 



(2.4) 
□ 



with 



Rh{p) 



-d 



I {\P?' - 1) dp ml + {2nh)-'' [ Ru{p) dp (2.5) 
J\p\<i -'IpKi 



{\p + + \p- - ) |0(7?)pd7? . 



We proceed to estimate Rh{p)- Note that for any a > 

max ((a + t)"" + (a - tf) = 20^ . 

\t\<a 

Taking a = |pp + |r/p and t = 2p ■ r] we deduce that 

1 (|p + ^|2s + |p _ ^|2s) _ |p|2s < (|^|2 ^ |^|2)s _ |p|2s ^ 

Next, for < s < 1 concavity implies that (a + b)^ < a* + sa^^^b for a, 6 > 0, from which 
we learn that 

(|p|2 + |^|2)s_|^|2s < ^|^|2(s-l)|^|2^ 

Hence, replacing ij with hr] and using (I1.14p we can estimate 

Rh{p) < s J \p\'^+^'\hr]f\^{r])fdr] = s\p\-^+^' j \V(l)\^dx < Ch'^\p\-^+'^' . 



Thus the upper bound follows from (|2.ip and (j2.5p . 
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3. ASYMPTOTICS ON THE HALF-SPACE 

Our goal in this section is to prove the analogue of Proposition [5] in the case where O is 
the half-space = {(x', Xi^) : > 0}. We define the operator on L^(R^J_), in the same 
way as Hfi, with form domain 

n'iRi) = |w G H'iR"^) : = on M'^ \ mJ| . 

We shall prove 

Proposition 8. Assume that cf) G Cq(M'^) is supported in a hall of radius I > and assume 
that is satisfied. Then for h > and any < 5i < 1 and < ^2 < min{l, 2s} we have 

< Tr {(l)H+(l))_ - L^^l [ (l)\x)dxh~^ + if I [ (l?{x', 0)dx'h-'^+^ 



a!-l-<5i T -a!+l+<5i 



This result depends on a more or less explicit diagonalization of the operator , which 
is far from obvious. This is accomplished in Subsections 13.11 and 13.21 relying crucially on 
recent results of Kwasnicki [KwalOa] about non-local operators on a half-line. These results 
are collected and extended to our needs in the appendix. 

3.1. The model operator on the half-line. In this subsection we collect some facts about 
the one-dimensional operator 



in L^(R+) with form domain and about the corresponding operator A in 

defined analogously to A'^ , but with form domain H^(R). 

For n > and t,u ^ M+, let e~^{t, u, //) and a^(i, u, fi) be the integral kernels of {A'^ — ^)^_ 
and {A'^ — f^)-, respectively. Similarly, we define a{t, u, fi) via (A — . To simplify notation 
we abbreviate a~^{t,fi) = a'^{t,t,fi). We also note that a(/i) = a{t,t,fj,) is independent of 
t G M+. The inequality A~^ > 1 implies that a'^{t,u,fj,) = e~^{t,u,fi) = for fi < 1 and 
similarly for a(t,u,iJ,) and e{t,u,n). 

The following two results about e'^{t,iJ,) and a~^{t,fi) are rather technical and we defer 
the proofs to Appendices IB. II and IB.2I The first one provides a rough a-priori bound on 
e+(t,M,//). 

Lemma 9. For any fi > and t,u G M_|_ one has |e^(t,u, < C^^^"^^ . 

The second result in this subsection quantifies that a~^{t,fi) is close to a{p) for large t. 
Lemma 10. For any < 7 < 1 there is a constant such that for all /U > 1, 

) 

\a+{t,n)-a{^)\dt<C^I^ {{In Ijf + 1) . (3.1) 
In particular, the function 

m = / \iT''{a{\i'\-^')-a+{t\aW\-''))di' , t>0, (3.2) 
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satisfies for every < 7 < 1 



poo 

/ t^\K{t)\dt < 00. 
Jo 



(2) 

With this lemma at hand we can now define the constant ^ which appears in our main 



theorem by 



= / K{t)dt. (3.3) 







(This integral converges by Lemma flOl) Expression (I3.3P suffices for the proof of our main 

-(2) 

^s4' 



result. In Section [6l see also (]B.5p . we will derive different representation for L i. 



3.2. Reduction from the half-space to the half-line. Our goal in this subsection is to 
write the spectral projections of the operator on the half-space in terms of those of the 
operator on the half-line. Before turning to spectral projections we treat resolvents. 

Lemma 11. For x = {x',X(i) G M^, y = {y',yd) G o-nd z £ C \ [0,oo) the resolvent 
kernels of and of A'^ are related by 

This lemma, together with the representations (see, e.g., |Kat66] ) 

^7r« Jr(/i~2») 

and 

iH+)- = -h^' / -z)-^dzdv, 

Jo 27ri Jy(u) 

where r(z/) = {z £ C : \z\ = 1^}, implies that 



and 

We now give the 

Proof of Lemma By scaling we may assume that h = l. Given / G L2(M^) we want to 
solve {-AYu = zu + f. Take V G C^{W^~^) and (j) G C^(0, 00). Then -0 (g) </> belongs to the 
form domain of and therefore the equation implies that 



{W? + edyHi')mm\id)di' did =zjj i^ie)m)H^',^d)dedid 

+ 11 4>{e)m)fie, Cd) ^e' dCd . 

Since ip is arbitrary, this means that for a.e. 

(le? + 6rM^uiC', dU = zj J{^)u{^', U) did + j J{U)f{C', U) did ■ (3.6) 
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For fixed ^' define functions v^i, g^i and XE,' o^i (0,oo) by 

v^,{t) := (27r)-('^-i)/2 I |^'|-iu(x', |e'r4)e-^«'-^' dx' , 



-1 



g^,{t) := (2vr)-(^-i)/2 / f{x' , \e\-'t)e-'^'-'' dx' , 

We note that the one-dimensional Fourier transforms of these functions are given by 

Hence we can rewrite (13.61) as 



J{l + kyxi'{k)v^'ik)dk = ler^^z J x^.{k)v^,{k)dk + I x^,[k)g^,{k)dk. 

Note that both V(^i and Xi' belong to the form domain of and that the set of all functions 
X^i obtained in this way is dense in the form sense (by the definition of A^). Therefore the 
equation can be written as 

We abbreviate rf{t,w) = (yl+ — z)~^{t,w) and conclude that 



v^'it) = J rj^^,\-2sit,w)g^f{w)dw . 
Recalling the definitions of v^' and g^' this reads 
\('\~'^u{x',\^'\-H)e-'^'-^' dx' 

rtu,-2s{t,w)\e\-'-''f{x',\e\-h)e-^^'-^' dwdx'. 



-1 

fOO 



Multiplying by |^'|, setting t = \S,'\xd and inverting the Fourier transform, we obtain 



= (27r)-'^+i / / r+ .,(|e'|xd,le'|yd)ICV"'V(y)e*^'■("'-^')^^e't^y• 
This proves the lemma. □ 

3.3. Proof of Proposition [8l Our next step is to state upper and lower bounds on 
Ti {(pH~^(j))_ in terms of the one-dimensional model operators A and , in particular, 
in terms of the function K{t) given in (|3.2|) . As explained below, the main result of this 
section, Proposition [51 will be a direct consequence of the following estimates. 
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Proposition 12. Assume that (j) G Cq(M°') is supported in a hall of radius I = 1 and assume 
that ^1.14\ ) is satisfied with 1 = 1. Then for any < 82 < min{l, 2s} there is a constant 
such that for all h > we have 



h \ h 



dxh 



(3.7) 



TV [4>H+(f)_ < L^^l [ c^'^{x)dxh~'^ - [ 

TV {(^H+(l))_ > L^p^ [ (t)\x)dxh-^ - I (f'{x)\K (^) dxh''^+^ - Cs,h~'^+^+^^ . (3.8) 
' Jm.<1 Jr'I h ^ h J 

Assuming Proposition 1121 we now give the short 

Proof of Proposition\^ To prove the proposition we may rescale (p and hence assume 1 = 1. 
Proposition [8] is then an immediate consequence of Proposition [12] provided we can show 
that for any < 5i < 1 there is a C^^ such that for all /i > 



^\x)\k T-^) dx - Lf\ j <p\x',0)dx' 



(3.9) 



In order to obtain the latter bound, we substitute Xd = th and write, recalling ()3.3p . 

th 



1 



4>\x)^K {^)dx- L''^^ I 4>\x', 0)dx' 
By Holder's inequality we can further estimate 

rth 



(2) 







drClr{x',T)dTdx'dt. 







R<*-i Jo 



drr{x',T)dTdx' 



< 



th \ "^i 
dT 







drrix', T)dx' 



dT 



< Ct^'h^' . 



Since t^'^\K{t)\ dt < 00 hy Lemma [TOI we obtain inequality ()3.9p . 



□ 



In the following two subsections we shall prove the lower and the upper bound in Propo- 
sition [T21 respectively. 

3.4. Lower bound. To prove (j3.7p we use that 

-Tr{(j)H+(j))_ > -Tr (0(i?+)_(/)) . 
The lower bound follows from this by integrating the identity 



{H+).{x,x) = h-^L^p-h-^K[^ 



(3.10) 



against cp"^. Equation (|3.1Up is a consequence of 
Subsection 13.21 we learn that 

{Ho)-{x,x) 



Indeed, by the same argument as in 



/|l+2s 



(27r)°'-i h'i 

On the other hand, by direct diagonalization as in Subsection 12.11 we find that 

iHoUx,x) = h-''L^p. 
Comparing these two identities with (|3.5p we arrive at (|3.10p . thus establishing (j3.7p . 
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3.5. Upper bound. To prove ()3.8p we set 7 = {H~^)^. Its integral kernel is given by (|3.4p 
in terms of the kernel e~^{-,-,fi) of {A~^ — By the variational principle it follows that 

-Ti {(j)H+(j))_ <Tr {(j)-f(j)H+) 
1 



h2d 




X (Ipl- - 1) e-(--)/'^ H^) m ^^^^j^ . (3.11) 
We insert the identity 



use the symmetry in x and y and substitute q = Pd/\p'\ to obtain 

-TV (077+0) < 4[0]-i?ft,[0] 

with the main term 



" —+ " — + " — " " 

and the remainder 

Jr^ Jr^ Jr-^-i JRd 

x|ppVH.-)M|0(,)-0(,)P^|lJ^. 
Since G Co(M'^) we can perform the ?/'-integration in //j We use the fact that 

p poo 

/ / e■^{xd,yd.^Ji){{q^ + lr -^l)e'^y^^-''^'^''dyddq=-a+{xd,Zd,^l) 
JR JO 

and obtain 

h'^^^ jRd Jo JRd-i Jr 



h"^ iRt Jr-^-i ^ ' IS , ; ^27r)rf-i ■ 

Using again p.lOp we find that 

h[cP] = -Lfl[ (l?[x)dxh-^+ f cp\x)K(^)dxh-^. (3.12) 

' JR'l JR| \ n / 

It remains to study Rh[<P]- We claim that for any ^ — s < a < min{|, 1 — s} there is a 
such that 

iRhM < ^313) 
for all /i > 0. This, together with (j3.12p will complete the proof of (|3.8p . 
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In order to show ()3.13p we perform the p integration and find that 



Rh 



C 



f f I leie 



h 



h 'le 



/|2s 



\x-y 



d+2s 



We insert 



jCix'-y')/h 



and integrate by parts to get 



jjd-2s—2(j 



l\l-2aii'ix'-y')/h + 



xd\i'\ ydW\ 



h '1^ 



l\2s 



\x-y 



\d+2s 



dx dy . 



By Lemma [9] and the fact that e~^{t, n, /x) =0 for fi < 1 we arrive at 



c 



< 



]^d-2s-2(T 

c 

J^d~2s-2cT 



:|5'I<1} 

\(t>{x) - (t>{y)\' 



(x)-</.(y)p 



\x -y 



d+2s 



dxdy 



\x-y 



d+2s 



dxdy . 



According to Lemma [24] this imphes ()3.13p and hence completes the proof of (13. Sh . 

4. Local asymptotics near the boundary 

In this section we prove Proposition [5j After having analyzed the half-space case in the 
previous section, we now show how the case of a general domain follows. We shall transform 



the operator locally to an operator given on the half-space 



{iy',yd) G 



yd > 0} and we shall quantify the error made by this straightening of the boundary. 

Under the conditions of Proposition \E\ let B denote the open ball of radius / > 0, con- 
taining the support of 4>. For xq € BD dQ let be the inner normal unit vector at xq. We 
choose a Cartesian coordinate system such that xq = and = (0, . . . , 0, 1), and we write 
X = {x', Xd) G M'^^i X M for X G W^. 

For sufficiently small I > one can introduce new local coordinates near the boundary. 
Let D denote the projection of B on the hyperplane given by x^ = 0. Since the boundary 
of Q is compact and C^'" there is a constant c > such that for < / < c we can find a real 
function / E C^'" given on D, satisfying 

dnnB = {{x',Xd) : x' e D,Xd = fix')} riB . 

The choice of coordinates implies /(O) = and V/(0) = 0. Hence, we can estimate 



sup \Vf{x') 
x'eD 



sup |V/(x')-V/(0)| < Cf\x- 
x'eD 



Since the boundary of is compact we can choose a constant C > 0, depending only on ri, 
in particular independent of /, such that the bound 



sup |V/(x')| < Cf 

x'&D 



(4.1) 



16 



RUPERT L. FRANK AND LEANDER GEISINGER 



holds. 

We introduce new local coordinates via the diffeomorphism ip : D x 



given by 



and 



Uj = ^Pj{x) = Xj for j = 1, . . . ,d — 1 



yd = ^Pd{x) = Xd- f{x') 



Note that the determinant of the Jacobian matrix of (p equals 1 and that the inverse of (/? is 
given on ran 9? = D x M. In particular, we get 

ip {dQ nB) c dR'i = {y eR"^ : yd = 0}. 

Fix V £ T-l^{Q) with support in B. For y £ ran 93 put v{y) = v o (p^^i^y) and extend v by 
zero to M"'. 

Lemma 13. The function v belongs to H'^{M.'i!^) and for < I < c we have 

^ we find Xd = yd+f{y') < f{x'), 



Proof. By definition, v belongs to H'^{W ) and for y G 
thus v{y) = v{x) = 0. Therefore v belongs to 'H^{M.'^ 
Using the new local coordinates we get 

\v{x) — v{w)\'^ 



{v,{-A)l,v) = Cs,d 



x — w 



\d+2s 



-dx dw = Cs d 



\x — w\ 



where y = (p{x) and z = (p{w), thus x = {y' ,yd + f{y')) and w = [z' ,Zd + f{z')). Let us write 
1 1 



1 



\y- z 



d+2s 



\y - z 



d+2s 



[W - z'? + {yd + f{y') - zd- f{z')?f^^' 

After multiplying out, the last fraction equals 

Uiy') - f{z')f + 2{yd - zd){f{y') - /(z')) \ "^'^'"''^ 



1 + 



\y - z\ 



and we can employ ()4.ip to estimate 

(/(y') - f{z')f + 2{yd - Zd)U{y') - fiz')) 



\y - 2| 

< sup I V/p 'f' ~ ^'f + 2 sup |V /I ~ ''^ ^'^l < Cf 



\y - z\ 



\y - z[ 



Choosing I small enough we can assume CI" < 1/2. Then, combining the foregoing relations, 
we find 
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< C,,, II \v{y)-v{z)\' 

< Cl"(v,(-A)Uv). 
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1 



1 



^|rf~|~"2s \x X/-^ I '^^^^^ 



dy dz 



This proves the first claim of the Lemma. The second claim follows by interchanging the 
roles of and (-A)f,. □ 



On the range of (p we define (f)u — </'u 

o if ^ and extend it by zero to M'^ such that 
lu S Cq(M'^) and ||V0u||oo < Cl^^ hold. Using Lemma[T3]we show the following relations. 



Lemma 14. For < I < c and any h > the estimate 



holds. Moreover, we have 



(j)^{x)dx = / <t)'^{y)dy 



and 



< / (l)\x)da{x) 
Jdn 



^\y',0)dy' < CI 



d-l+2a 



(4.4) 
(4.5) 
(4.6) 



Proof. The definition of (j) and the fact that the Jacobian of (f) equals 1 immediately gives 
3|) . Using (j4.ip we estimate 



0^(x) da{x) 



02(y',O)Vl + |V/|2dy' < / 4>\y',o)dy' + ci^~'+'". 



an 

from which ()4.6p follows. 

To prove (j4.4p we refer to the variational principle once more and note that 

-Tj:{(l)Hn<P)^ = mi Tv {<p-f<pHn) , 

0<7<1 

where we can assume that infimum is taken over trial density matrices 7 supported in B x B. 
Fix such a 7. For y and z from D xM set 

i{y,z) = 7 {'P^\y),'P'Hz)) , 

so that < 7 < 1 and the range of 7 belongs to the form domain of (pH^cj). According to 
Lemma [T3] it follows that 

Tr (<A7'A^f7) > TV U^^ (h^'H - Cr)(-A)^, - 1 



and consequently 



Tr (0Fn0)_ < Tr ( ( (1 - Cnh^'{-A)', - 1 
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Set e = 2Cl" and assume / to be sufficiently small, so that < e < 1/2. Then 
'TT{cPHn^)_ < Tr(<^((l-Cr)/i2«(-A)^, 

< Tr ((-/i'A)^, _ i) 0) + Tr ((e - C/")/i2«(- A)^, " ^) „ 

< TV(<^i7+<^)_ + eTr ((/i272)(-A)^, " ^) _ ■ 

Using Lemma El we estimate Tr {4){{h^y 2) {-Ay^^ - l)"^)- < Cl'^h~'^ and it follows that 

Finally, by interchanging the roles of and H^, we get an analogous lower bound and the 
proof of the Lemma is complete. □ 

We conclude this section by giving the short 

Proof of Proposition 0. It suffices to combine Lemma [H] and Proposition [H □ 

5. Localization 

In this section we construct the family of localization functions {(pu)ueR'i prove Propo- 
sition [3l Fix a real-valued function cf) G C^(R'^) with support in the ball {x G M'^ : |x| < 1} 
that satisfies \\(p\\2 = 1- We recall the definition of the local length scale l{u) from (jl.lip . 
For n, X E M"^ let J{x, u) be the Jacobian of the map n 1— )• (x — u)/l{u). We define 

such that 4)u is supported in the ball Bu = {x : \x — u\ < Hu)}. 

By definition, the function l{u) is smooth and satisfies < l{u) < 1/2 and ||V/||j^ < 1/2. 
Therefore, according to jSS03| . the functions (pu satisfy (|1.12p and (|1.13p for all u G M.'^. 

To prove the lower bound in Proposition [3] we follow some ideas from |LY88J. In particular, 
we need the following auxiliary results; the first one gives an IMS-type localization formula 
for the fractional Laplacian. 

Lemma 15. For the family of functions (0u)^g]gd introduced above and for all f G ^{^{O,) 
the identity 

(/,(-Ar/)= / icPuf,{-ArKf)iiur''du-{f,Lf) 

holds with $7* = {n G : suppc^u n 17 7^ 0}. The operator L is of the form 

L = f L^Jiuydu, (5.1) 
Jn* 

where L^^ is a bounded operator with integral kernel 
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Here denotes the characteristic function of il. 
Lemma [T5] imphes that for any operator 7 with range in 

Tr7(-A)^ = / TT{-i(t)u{-^y<Pu)l{u)-'^du-TT^L. (5.2) 

The next result ahows to estimate the locahzation error Tr7L. 
Lemma 16. For u € M"^ and < 5 < 1/2 we have 

Tr7L<^„ < TT^{C5^-^%u)-^'xsXn)+CU\l{u)-^'5-''+^-^'T{5) 

with 

1 if 1 - d/4 < s < 1 

r{5) = I |hi(5| if < s = l-d/4 . 

jd+4s-4 if < s < 1 - (i/4 

where xs denotes the characteristic function of {x G M"^ : \x — u\ < l{u){l + (5)}. 

Proof. By translation and scaling we can assume that u — and ^(w) — 1, and hence (p^ — 0* 
(This rescaling changes but the bound we are going to prove is independent of the domain 
and therefore not affected by this dilation.) We set 



L,j>{x,y)x5{x)x5{y) ii\x-y\<5 
if |a; — y| > (5 



L^{x,y) = L^{x,y) — L^{x,y) and 6{x) = f L^{x,y)dy. By a simple adaption of the 
arguments of |LY881 Thm. 10] we find that for any e > 

Tr7L^ < Tr7(0 + exo) + ^Tr(LO)' . (5.3) 

It remains to bound 6 and Tr(L^)^. 

We begin by estimating 9. By definition, for |x| > 1 + 5 we have L^{x,y) = and hence 
9{x) = 0, and for \x\ <1 + 5 we get 

0{x) = a,[ ^ ^-^}^^^^dy<C\mt [ , IjT^dy. 

\y\<l+5 ' 

Thus, for all x£W^ 

9{x) < C 5^-^' Xs{x) . (5.4) 
Finally, we estimate Tr(L^)^. The symmetry of L^{x,y) implies 

^2x 2 



where A denotes the set {{x,y) € M'^ x M*^ : |x| < min(|y|,l), |x — y| > 6}. Set Ai = 
{{x,y) e A : \y\>2} and A2 = {{x,y) e A : \y\ < 2}. Then 

Tr(L°/ < C I [^^ ( ^^ ^y^^^^ ^ dxdy + C J |^ _ ^,1+4.-4 dxdy 

< Cr{6). (5.5) 
Choosing e = 5^^^'^ and combining (|5.3p with (|5.4p and (j5.5p yields the claimed result. □ 
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Proof of Proposition\^ We apply Lemma [16] with a parameter < 5^ < 1/2 to be specified 
later. For ease of notation we write Xu instead of xSu- Identities ()5.ip and (j5.2p and the 
estimate from Lemma [TBI imply 



Tr7(-Ar > / TT^{<l^u{-^)^u-C5l-^n{u)-^\u{x)xn{x))l{ur''du 

1* 

ChW [ 5-'^+^'^'r{5u)l{u)-'^~^' du. (5.6) 



If the supports of Xu and overlap, we have \u — u'\ < {3/2)l{u) + l{u'). It follows 
that l{u') - l{u) < \\Vl\\oo {{3/2)l{u) + l{u')). Since ||VZ||oo < 1/2 we find l{u') < Cl{u) and 
l{u)~^ < Cl{u')~^. Similarly, we get l{u) < Cl{u'). We assume now that 6^ satisfies 

Su<C6u' ii\u-u'\<{3/2){l{u)+l{u')). (5.7) 

Using these locally uniform bounds on l[u)/l{u') and 6^/^^', together with (jl.l3p . we can 
deduce the pointwise bound for all x S M"^ 



Sl-'niu)-'^Xu{x)xnix) j^ = l^^ 6l~'%ur'^xu{x)xn{x) {^j <t>l{x)-^^ 



du 



2-2su..'\-2s ^ du' 

l{u'Y ' 



<C / 4>u'{x)dt7'n{u'y<p^,{x) 
Jn* 



Rewriting the last integral with u as integration variable, in view of (15. 6p . we find 



By the variational principle it follows that 

TV(lfn)- = -^mf Jr7((-/i2A)^-l) 



Jn* 



-2s\ , \ du 



l(uY 



du 



+C.-y^^5-.«-V(^.)^. (5.8) 

To bound the first term, we use Lemma [6l For any u G M*^, let pu be another parameter 
satisfying < pu < 1/2 and estimate 

Tr (0„ [{-h^/\y - 1 - Ch^'5l~^' l{u)-^') <Pu)_ 

< Tr {cl)uHn<Pu)^ + CTr (0„ {p^h^'{-Ar - - h"' 5l'^n{u)~^') 

< Tr {ct>uHM^ + ci{ur{puh'T'/^'^^ {pu + h^'sl-'niuy^f^"/^''^ . 

We pick pu = h?^ d"^"^^ [{u)""^^ . By (|1.17p and our assumption that 5u < 1/2, we see that 
Pu < (/i//o)^''2^*~^. We assume now that h < C~^lo (with a possibly large constant C) in 
order to guarantee that pu < 1/2. With this choice we find 

Tr U J {-h'AY - 1 - „ " U J <Tr [cj^uHnM- + C " J_J • (5.9) 
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Combining ()5.8p and ()5.9p we obtain 

du 



Tv{Hn)- < 



+ C 



+ l^n- (5.10) 



At this point we choose 6u in order to minimize the second integrand, which we shah 
denote by /„. We pick 



h/l{u) if 1 - d/4 < s < 1 

{h/l{u))\ ln(/(n)//i) 1 1/(4-4^) if < s = 1 - d/4 
{h/l{u)fl^*-^'^ if0<s<l-(i/4 



and note that < 1/2 if /i < C ^/q by (|1.17p . Moreover, (|5.7p is an easy consequence of 
the corresponding estimate for l{u)/l{u'). With this choice we arrive at the bounds 



Iu<C{ 



h-'^+^l{u)-'^ if 1 - d/4 < s < 1 

/i-'^+2/(u)-2|ln(/(n)//i)|i/2 if0<s = l-d/4 • 

^-d/2+2si^^yd/2^2s if < S < 1 - (i/4 

Finally, we integrate with respect to u. The same arguments that lead to (|1.18p and (|1.19p 
yield 

^ /i-'^+^Z-i if 1 - d/4 < s < 1 

/i-'^+2;-i|ln(/o//i)|^/2 if0<s = l-(i/4 . 

h-d/2+2si^d/2-2s+l if < s < 1 - (i/4 

This completes the proof of the lower bound with the remainder stated in Proposition [3j 
To prove the upper bound we put 

4)u {4>uHn(puf_ (puKuy^ du. 



ludu < C < 



Obviously, 7 > holds and in view of ()1.13p also 7 < 1. The range of 7 belongs to 
and by the variational principle it follows that 

-Tr(Fn)- <TTiHn = - I Tr /(n)"'^ dtx . 

This yields the upper bound and finishes the proof of Proposition [3l 



□ 



6. Discussion of the second term 



6.1. Representations for the second constant. In this section we study the second term 
of (|1.6p in more detail. First we derive representation (jl.Sp . 



Proposition 17. One has 



L 



(2) 
s,d 



d-l 

;d-2| 



aw 



/|-2s 



dp' 



(27r)^-i 
2s 



(27r)'^-i {d - l){d - I + 2s) 



Tr 



-{d-l)/2s 



X-{A 



+ \-{d-l)/2s 



(6.1) 
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Here x is the characteristic function of M-|- and 

/•oo 

C(/")=^"W {aifi)-a+it,fi)) dt. (6.2) 
Jo 

Proof. The first identity follows immediately from (I3.2p and (j3.3p . The second identity 
follows from the fact that 

y\^'(E -\v'\-^') ^ j^t^ ?f rj-id-i)/2s 

i-i'-^' ^ l-^' ^"(27r)'^-i (2^)°'-i (d-l)(d-l + 2s) 
for any E > 0, which by the spectral theorem implies that 

and similarly for A. □ 
Remark 1. There is another representation, namely, 

where 

^(;u)= / (e(/.)-e+(t,/.)) dt. (6.4) 

JO 

Here e(//) and e^{t,fi) are the diagonals of the integral kernels of the spectral projectors 
{A—n)^ and {A~^ — fj,)^, respectively. We have not shown that the integral in (j6.4p converges, 
since we will not use (|6.3p in the remainder of this paper. Identity (|6.3p is an easy consequence 
of ()6.ip and the fact that 

o-ifA — I e{T) dr a'^{t,fi)= / e'^{t,T)dT 
Jo Jo 

which follows by the spectral theorem from {E — ^)_ = J^{E — t)^ dr. Representation (16.30 

is natural since in terms of this function the conjectured formula for the number of negative 

eigenvalues of takes the form 

2, ^so ^pdx ff ,fi„n-2s.dp^da{x) ^^^^ 



which is the analogue of well-known two-term semi-classical formulas in the local case; see, 
for instance, |Ivr801[S^V96j . The function ^ plays the role of a spectral shift. Note that we 
avoided to write (|6.2p and (j6.4p in terms of a trace. While the integrals on the diagonals 
converge, we do not expect the operators to be trace class, see |Pu08) . 

Remark 2. Yet another representation is 

/■oo 

^"'+'41 = / {H-{^,x) - (i7+)_(x,x)) dxd. 



(Note that the right side is independent of x'.) This follows from (j3.5p and the corresponding 
formula for H. Using this representation one sees that our asymptotic formula coincides with 
the one obtained in IBK08l[BK509l . 



Finally, we refer to ()B.5P in the appendix for a representation of L], ^ in terms of generalized 
eigenfunctions of A'^ . 
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6.2. Positivity of the constant. Here we shall prove 



(2) 



Proposition 18. For any < s < 1 and d > 2, one has ^ > 0. 
Proof. We use the second representation in (|6.ip and the fact that 

for every > to see that 

lsrf-2| 



T{{d-l)/2s) 



-l3Ep{d-l)/(2s)-l^P 







L 



(2) 
s.d 



2s 



(2^)°'-i {d - l){d - 1 + 2s)r((d - l)/2s) 



oo roo 



JO 



By means of the Trotter's product formula (see also jBK08j for a probabilistic derivation) it 
is easy to see that 

exp{-f3A)(t,u) > exp(-/3A+)(t,n) 
for every /3 > and every t,u > 0. This, together with the fact that A ^ A'^, proves the 
proposition. □ 

6.3. Comparison with a fractional power of the Dirichlet Laplacian. It is well- 
known that the Dirichlet Laplacian — Aq on satisfies 

Tr {-h^An - 1)_ = \n\ h^'' - if^ \dn\ /i-'^+i + o(/i-^+i) , 

see, e.g., [FGllj for a proof under the sole assumption that dil G C^'" for some < a < 1. 
Here 



L 



(1) 



^ / {\P\'-I).dp 



i'"^ " (27r) 

and, by an argument similar to that in our Proposition 1 171 one can bring the second constant 
in the form 



L 



(2) 



2d-2\ 



(27r)'^-i {d-l){d+l) 



where B 



d^ 



+ 1 in L2 



and 



Tr 



d^ 

dp 



XB 



-{d-l)/2 



-{d-l)/2 



+ 1 with Dirichlet boundary conditions in 



L (M-|_). A short computation, using the fact that 



{E' - 1)_ = s(l - s) [\e - t)-t'-^ dT + s{E - 1)_ , 
Jo 



gives 



- lS'] \dn\ h-^+^ s (^(1 - s) r('^-i)/2+-i dT + ?j+ o{h-^+') 



that is, 



^d-2l 



d-l + 2s 
2s 



~A2) ^ S{d+l) (2) ^ 

'^'^ d-l + 2s i''^ (27r)'^-i (d- l + 2s) 



Tr 



xB 



<d~l)/2 



X-{B' 



-{d-l)/2 
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Since 



we find that 



im _ (2) ^ IS'^-I 2. 

s,d (27r)rf-i (d- l + 2s) 



We now apply Lemma [T9l below with B = —(f'/dt^ + 1 in with P being the projector 

onto L2(M+) and with (t){E) = E". Then (f){PBP) = {B+Y and P(I){B)P = yl+, and 
therefore (|6.5p yields 

(5+)" > A+ . 

Since E i— t- E^^'^^^^^'^^ is strictly monotone and since the operators A'^ and {B'^y are not 
identical, we conclude that 



Tr ((^+)-(^-i)/2. _ (^+)-(d-i)/2^ > o_ 



This shows that ^ — ^ > and completes the proof of Proposition [2j □ 

In the previous proof we used 

Lemma 19. Let B be a non-negative operator with ker B = {0} and let P be an orthogonal 
projection. Then for any complete Bernstein function 

(piPBP) > P(j){B)P . (6.5) 

We recall (see, e.g., [SSV10| ) that complete Bernstein functions (also known as operator- 
monotone functions) are characterized by the representation 

poo 

with a G M, 5 > and a positive measure p satisfying 



-dp(r) (6.6) 



Jo ^(TT7)^^(^)<°°- 



Proof. We first prove that 

B-^ > {PBPy^ . (6.7) 

Let us write P^ = 1 — P, so that 

B = PBP + P^BP^ + P^BP + PBP^ . 
By the Schwarz inequality we have 



for any e > 0, and hence 



P^BP + PBP^ < ePBP + -P^BP^ 



B<{1 + e)PBP + ( 1 + ^ ) P^BP^ 



Since P and P^ are orthogonal we can invert this inequality and obtain 

B-^ > -^(PBP)-^ + -^(P^BP^)-^ . 
Thus dlZl) follows by taking e 0+. 
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Now if (j) is of the form (16.60 then 

POO 

(p{PBP) - P(j){B)P = - {{PBP + r)-^ - P{B + tY^'P) dp{T) . 

Jo 

By (16. 7p with B replaced by -B + r, this is non-negative. □ 



Appendix A. Equivalence of (II. 2p and (|1.3p 
For the sake of completeness we include a short proof of 

Lemma 20. Let {Xk)km be a non- decreasing sequence of real numbers and let A,C > 0, 
B,D gM. and —1 <a — l<b<a be related by 

C = .4-i/"a(a + , D = B{A{a + . 

Then the asymptotic formula 

N 

Y,>^k = AN^+^+BN''^\l + o{l)), N^oo, (A.l) 

k=l 

is equivalent to 

^(A-Afc)+ = CA(^+°)/''-Z)A(^+*)/°(1 + o(1)), A^oo. (A.2) 

ken 

Proof. This lemma is a consequence of Hardy, Littlewood and Polya's majorization theorem, 
which says that for any non-decreasing sequences {ofc} and {6^} 

N N 

^ Ofc < ^ 6fc for all G N (A.3) 

k=l k=l 

is equivalent to 

oo oo 

^(A - afc)+ < ^(A - bk)+ for aU A G M ; 

fe=l k=l 

see, e.g., |M079l Prop. 4.B.4]. As usual, we will denote property (|A.3j) by {ak} ^ {bk}- 

We fix e > and set = A{a + l)k"' + {B ± e){b + 1)A;^ Note that the assumptions on 
a and b imply 

N 

Y,Pt = + {B± e)N''+\l + 0(1)) , A ^ oo , (A.4) 

k=l 

and 



(A(a + l))i+Va (A(a + l))(i+'')/'^ 



(A.5) 



First, we assume that ([AT]) holds. Then, by (IATT]) and (IXil) there is an Ae G N such that 
for all N>N^ 

N N N 

fc=l k=l k=l 
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We put = (3^ for k > and = max(/3^, A^), = min(/3^, Afc) for k < N^. Thus 
{a^} -< {Afc} -< {Ofc }, and therefore 

^(A - «+)+ < ^(A - Afc)+ < j;(A - «-)+ for all A G R . 

fceN fcGN keN 

Since ZlfceNl^ " "fc )+ = 'Ekeni^ " + the assertfon ([A^]) follows from (1X5]) . 

The converse implication is proved similarly. □ 

Appendix B. The one-dimensional model operator 

Here we outline the calculations that are necessary to complete the analysis of the model 
operator introduced in Section [3l The results depend on the following spectral represen- 
tation of the operator found in [KwalOa] . 

Theorem 21. For E > let 

tpiE) = {E + iy-i 

and for A > put 7a(0 = Oi/0<^<l and 

1 Xij'iX"^) sm{TTs){C'^ -ly 



vr V^(A2)2 + (^2 _ ly _ 2^(A2)(e2 _ i) cos(7rs) 



if i> 1. Moreover, define a phase- shift 



and functions 

l-OO 

= sin(Ax + i?A)+ / e-"«7A(0de, x>0. (B.2) 

JO 



Then 



oo 



^f{\) = \l-l f{x)Fxix)dx 



IT 

defines a unitary operator from L2(M_|_) to L'^ 

This operator diagonalizes in the sense that a function f £ L2(M_|_) is in the domain 
of if and only if (A2 -|- 1)''<I>/(A) is in L'^iJR^), and in this case 

<^>A+f{x) = (a2 + ly-^fix) . 

According to [KwalOa] the Laplace transform of 7a is a completely monotone function 
bounded by one. Prom (|B.2p it follows that for all t > 

\Fx{t)\<2. (B.3) 

Theorem [21] states that the functions Fx are generalized eigenfunctions of the operator 
A~^. Hence, we can write 

e+(i,n,/x) = - / ((A2 + l)^-^)VA(t)FA(n)dA. (B.4) 
Jo 
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From ([33]) and (fL8|) it follows that 



i (l-2F,^(*))(A= + l)-<->/'.A.*. ,B,5) 

B.l. Proof of Lemma[9l Lemma[9]is an immediate consequence of ()B.4p . In view of ()B.3P 
we estimate 

\e+{t,u,fi)\ < C dX < (:7//(2s)_ 

This proves the lemma. 

B.2. Proof of Lemma IIOL First we need the following technical result about 'dx. 

Lemma 22. The phase-shift "dx is monotone increasing and twice differentiable in X > 0. 
It satisfies 

■do = and i^a — >■ ^(1 — s) as A — )• oo . 
The first and second derivatives are bounded and one has, as A — )• oo, 

h dHx 



d^x d^^x ^ ( 1 



dX dX^ ^ \X 

Proof. Following [KwalOaj . we substitute ( = Xz ior ( £ (0, 1) and C = A/z for C £ (liC«) 
in the definition of 'dx and obtain 



IT Jo 1 — \Z 

Note that the function 

1 V'(A2) - V'(A^-z^) 1 {l + X"^)' -{1 + X'^z^y 



iPixyz^) - v(A2) z^ (1 + xyz^y - (i + a2)- 

equals 1 for A = and that for all z G (0, 1) it is increasing in A > and tends to as 
A tends to infinity. By Lebesgue's dominated convergence we find 'do = and 

hm = I rr^Hz^'-^)dz = 



By ()B.ip . we also have 



1 1"^ 
^x = - bxiOdC 
71" Jo 



with 

A f sji+x^y-^jx^-e y 

C2-A2 V (A2 + 1)- - ((2 + 1) 

We remark that 

1 ,^ f se 

[i-ey-i, 

for all C G (0, oo). Since the last expression is integrable in € (0, oo) it follows that 



\dxbx{0\ < dxbx{C)\x=o = ^^^(^ 



= - dxbxiQdX 

dX TT Jo 



28 RUPERT L. FRANK AND LEANDER GEISINGER 



is bounded and, in particular, we obtain 

dX 

Similarly, we can show existence and boundedness of the second derivative and decay of the 
derivatives as A — )• oo by explicit calculations and Lebesgue's dominated convergence. □ 

To simplify notation we put 

, 1-^/A2 

- T^mNW) 

for E > 0. Moreover, we write G\ for the Laplace transform of 7a and g\ for the Laplace 
transform of Gx. According to [KwalOa] we have 



_ Acos^A+tsin^A .2 #(A2) Vxjt) 



with 



Mt) = exp ^- In (Va(C')) dC 

To prove Lemma [10] we need the following properties of (/3a- 



Lemma 23. The function 1 1— )• ipx{t) is differentiahle in t > and its derivative satisfies 

^'xi^) — ^i^) as A — oo , 



+ 0(A) asA^O. 

A=0 



Proof. For fixed C £ (0,oo) the function A i— )• ipxiC'^) is non-increasing in A > and tends to 
1 as A — )• oo. Moreover, 

C2 ^^"^^ V(c^ + i)'-iy 

is integrable with respect to ^ G (0,oo). Hence we find that 

1 r°° 1 

V'xm = - ^ln{i;x{e))dC 
^ Jo '^^ 

and f'xi^) = o(l) as A —7- oo by Lebesgue's theorem. 
In view of (|B.6p 

The second claim now follows from the fact that the derivative of A i— )• ^'x(^) bounded. □ 
Proof of Lemma \1(A In view of Theorem [21] we can write 

a(/i)-a+(f,/.) = - / {{X' + ir-f^)_{l-2Fi{t))dX 
^ Jo 

and by 

1 - 2Fx{tf = cos(2Ai + 2i?a) - 4sin(Ai + i^a) Gx{t) - 2Gx{tf . 



A=0 
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We get 



f'\a{iM) - a+{t,n)\dt < Ri{ix) + R2{n) 



with 



{/.v»_i)y2 



(/i - (A^ + 1)') cos(2At + 2t?A) dX 



dt, 



(fi - (A2 + ly) (2 sm(At + ^x) Gxit) + Gxitf) dX 



dt. 



To estimate Ri{^) we split the integration in t and integrate over t G [0,1] first. We 
assume < 7 < 1. The proof for 7 = follows similarly. 
We write 

Id , ^ „ , cos(2At + 2t?A) d^x 

and insert this identity in the expression for After integrating by parts in the A- 

integral one can estimate 

i^i - (A^ + 1)') cos(2At + 2t?a) dX 



dt < C^i ((In + 1) . 



To estimate the integral over t E [l,oo) we proceed similarly. We integrate by parts twice 
and get 

{fi - (A^ + 1)') cos(2At + 2'&x) dX dt < Cfi{lnfi + 1) . 



We conclude 

Riil^) < C^((ln/i)2 + l) 

and turn to estimating R2{^j)- 

Since Gx is non-negative and uniformly bounded, we have 

R2{fJ-) < C / (^-(A^ + l)^) / t"/ Gx{t)dudX. 

Jo Jo 

Identity ^Bl}i implies Gxit)dt = gx{0) and f^tGx{t)dt = g'^{0). We note that 

cos 19 X 



(B. 



5a(0) 



A 



Ijj'jX-^) 

V^(A2) 



and apply Lemma [22l to estimate Gx{t)dt < C (A A A ^). Moreover, by ()B.7p . 



5^(0) 



sini?A /V''(A2) , 



A2 



V'(A2) 



V'l(O) 



and we apply Lemma [22] and Lemma [23] to estimate tGx{t)dt < C (l A A It follows 
that 



/ t'^Gx{t)dt < c[l^x-^] 
Jo 
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Thus, by (IB.Sp . we arrive at 



i?2(^) < C 



(/i - (A^ + 1)*) (1 A A~^) dX < C^i (In /X + 1) . 



This finishes the first part of the proof of Lemma [101 
In order to prove the assertion about K{t), we bound 

/ t^\K{t)\dt< leT+^M t^\a+{t\£,'\,\£,'\-^')-a{\£,'\-^')\dtdi' . 

Jo •^I5'l<i -^0 

Here we also used that, since = a~^{t,^) = for < 1, we can restrict the integration 
in the definition of X to < 1. On the other hand, from (|3.ip we know that 

t^\a+{tfi-^/^', /i) - a(/x)| dt < C^/ii+(^+i)/(2^) ((In fif + l) . 

Combining these two bounds and using that 7<l<d — Iwe obtain the second part of 
Lemma [TOj □ 

B.3. A remainder estimate. The following technical lemma was needed in the proof of 
the upper bound near the boundary. 

Lemma 24. Assume that cf) E Cq(M'^) is supported in a ball of radius I = 1 and that (TTT^ 
is satisfied with 1 = 1. Then for any ^ — s < a < min{^, 1 — s} one has 



(-A, 



\x-y 



\d+2s 



dxdy < C 



(B.9) 



Proof. For x = (x', xa) G K'^^^ x E and y = (y', ya) € M"^"^ x E put 



{(l){x',xd) - (t>{y',yd)Y 



{\x'-y'\'' + {xd-yd?f^^' 
To establish (lR9ll we use that 

\m-<i>{y)? 



\x-y 



d+2s 



< c 



Fxd,y{x ) ^xd,y{^ )l 1 I 

\x' — 2'M-l+2o- 



(B.IO) 



and split the integration in x G M'^ and y G M'^ in four parts. First we assume that x and y 
are in Bi. Then we have to show that 

\Fxd,y{x ) ~ Fxa,y{^ ) 



\x' — _2'|'^-l+2o 



■ dz' dx dy 



Bi JBi J\x'~z'\<\x~y\/2 



\^--^d,y{x"^ Pxa,y{z')\ , 

\x' - ^'1^-1+2(7 ^ 

\^--i^y^'^')-^-'i^v^^')\.dz'dxdy 



lBiJBiJ\x'~z'\>\x-y\/2 \x' - z']"^ 1+2^^ 

is bounded from above. 

To estimate the first integral over |x' — z'| < |x — y|/2 we use that 

Hz') - Fix') = g 



(B.ll) 



i=i 



X' — Z' 



dt. 
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For j = 1 , . . . , d — 1 we have 



thus 

Hence, we obtain 



|(9,F.„,)(xO|<C|x-y|-'^+^-2^ 



|x - y\ 



< C\x' - z' 



x' + t 



[z' — x') 



2 Y 

+ {xd - Vdf dt 



l-a 



(B.12) 



with < a < 1 and /3 = {— — h s)/(a — 1), by applying Holder's inequality. Note that 
, , (z' — x') , x9 , ,9 9 fx' — -u') • (-z' — 3;') 

x'-y' + tY-, K +{xd-ydf = \x-y\^ +1" + 2t^ rr^^ 

\x' — z'\ \x' — z'\ 

> i\x-y\-tf . 

Inserting this into ()B.12p we get for |x' — z'l < \x — y\/2 

- F.,.,(xO| <C\x'- zT U {\x-y\- tf^dt 

<C\x' -zT\x-y\^^^+^^^^-''\ 
where (2/3 + 1)(1 — a) = — d — 2s + 2 — a. We conclude that for any 2a < a < \ and a < 1 — s 



\Fx^,y{x') Fx^^y{z') 

I I I, , I™/ -,/\d-l+2a 

' Bi J Bi J\x'-z'\<\x-y\/2 \X — Z \ 



■ dz' dx dy 



<C \x' -z'\~''+^-^"+''dz'\x-y\-''-^'+^~''dxdy 

JBi JBi J\x''^z'\<\x-y\/2 

<C. (B.13) 
Now we turn to the second integral in (jB.lip over |x' — z'| >\x — y\/2. Since 

0<F,„,(x')<|x-y|-'^-2^+2 (B.14) 

and o" < 1 — s we have 



^Xd, y{x') 
'Bi JBi J\x'-z'\>\x-y\/2 \x' — Z 

Moreover, 



'^^dz'dxdy<C [ [ 

I J B\ J B\ 



X — y|rf+2s-2+2fT 



< C. (B.15) 



J\x'-z'\>\x-y\/2 \X - Z \ 



1/g 



\x'-z'\>\x-y\/2 
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with ^ + ^ = li by Holder's inequality. Since a > ^ — s we can choose p > and 



p <? 

? > d^^- By (iBlil) . we have 



\x'-z'\>\x~y\/2 



It follows that 



\ 1/9 

Fl (z')dz' < C 



( f - VT + {Xd - Vdf) 



2^-qid/2+s-l) 



1/g 



<C\xd- yd\ 



-d-2s+2+((i-l)/ij 



Bi JBi J|x'-z'|>|x-j/|/2 



dz'dxdy 

I _ 7/ d— 1 + 2(7 ^ 



\X' — Z 



< C 



\X - y\-d+l-2aHd-l)/p _ y^\~d~2s+2+(d-l)/q 



< C t 



Bi JBi 
2 



d-2s+2+{d~l)/q I ^d-2 1-^2 ^ ^2N(-'^+l-2'7)/2+(d-l)/(2p) 



dr dt , 



where we substituted t = \xd — yd\ and r = \x' — y'\. Since p > and a < 1 — s we find 
f f f _J}^yi^') 

Jbi Jbi J\x'-z'\>\x-y\/2 W — Z' 

The estimates ([RTS]) and ([BJ6|) show that 

\Fxa,,y{x'') — Fx^^y{z' 



^;|i!i+2a dz'dxdy < C j\^-^'~^''dt < C. (B.16) 



Bi JBi J\x'-z'\>\x-y\/2 F ~ ^ I 



■ dz' dxdy < C 



(B.17) 



and from (jRTOl) . (|Bl3|) . and (iBTfl) it follows that 

,,|c/>(:E)-<^(y)|2 



\x-y 



\d+2s 



dxdy < C . 



The proof that the respective integrals over Bi x (M" \ -Bi), (M" \ Bi) x Bi, and (M \ 
Bi) X are finite is similar but easier, since supp0 C Bi and we only have to handle 

one singularity at a time. □ 
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